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Scope and Technical Challenge 

> Predict the strength and service life of ceramic & composite structures 

CMC - Ceramic Matrix Composites & PMC - Polymer Matrix Composites 

• Need to account for: 

- Wide variability in the strength of individual components 
(probabilistic/stochastic strength) 

- How strength changes with different types of loading 
(strength vs: multiaxial loading) and size of the structure 
(size-effect) 

- How strength degrades with time and fluctuating load 

- How strength/damage response of monolithic, anisotropic 
and composite material (architectures) differ 
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Approach /Outline 

1 . Overview: Describe the MAC and CARES codes 

- MAC/GMC: composite micromechanics model 

- CARES Unit Sphere: multiaxial stochastic strength model 

(isotropy & anisotropy ) 

2. Applying CARES to the MAC code to simulate stochastic damage 
progression in a ceramic matrix composite (CMC) 

❖ Cellular Automaton: Encouraging failure of adjacent elements - mimics crack-like growth 

❖ Visualization of element-by-element failure propagation for fiber, matrix, and interface 

- Status & Capability: Current progress of code integration effort 

> Examples: 

(1) Stress-strain response of a SiC-RBSN laminate (circa 1990) 

(2) Time-dependent degradation notional example 
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MAC/GMC Methodology: Generalized Method of Cells ( GMC) & 
High-Fidelity Generalized Method of Cells ( HFGMC) 


❖ 


Mjcimiechanjcs links the size scales & provides the composite response based on 


the composite constituent materials 
FEAMAC: MAC/GMC embedded in FEA as constitutive material 


Repeating Unit Cell (RUC) 
of composite material 


HFGMC (2000s) 

• 2 nd order displacement 
field in subcells 

• Elastic stresses and 
strains piecewise linear 

• Number of linear algebraic 
equations is rather large 

• Local inelasticity/damage 

• Has shear coupling 

• Has “subcell mesh ” 
sensitivity 

We currently only use GMC in FEAMAC/CARES 


GMC (1990s) 

• 1 st order displacement field in 
subcells 

• Stresses and strains piecewise 
constant 

• Number of linear algebraic 
equations function of number 
of subcells 

• Local inelasticity/damage 

• No shear coupling 

• No “subcell mesh ” sensitivity 


❖ RUC made subcells 

❖ Multiscale capability 



Aboudi, J.; Arnold, S.M.; and Bednarcyk, B.A. (2013) 
Micromechanics of Composite Materials: A Generalized Multiscale 
Analysis Approach, Elsevier, Oxford, UK. 


Aboudi, J; Pindera, M.J.; and Arnold, S.M. (2003): Higher-Order 
Theory for Periodic Multiphase Materials With Inelastic Phases. Int. 
J. Plast., vol. 19, pp. 805-847. 



CARES: Ceramics Analysis and Reliability Evaluation of Structures 


Life Prediction & Component Design Code For 
Advanced Ceramics 

• Developed to predict the probability of failure of ceramic 
components under complex thermomechanical loading 

• Combines Weibull & Weakest Link theory with 
concepts of linear elastic fracture mechanics (the 
Unit Sphere model) 



(CARES) 

reliability analysis 


Structural 

Model 


Element 

integration 

Point 


Component Reliability Analysis Capability: 

> Transient loads and temperatures 

> Fast-Fracture Rupture 

> Time-dependent (da/dt) crack growth 

> Cycle-dependent (da/dn) crack growth 

> Multiaxial stress failure models 

(PI A & Unit Sphere & Tsai-Wu & Tsai-Hill) 

> Proof test 
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Nemeth, Jadaan, Gyekenyesi.: “Lifetime Reliability 
Prediction of Ceramic Structures Under Transient 
Thermomechanical Loads.” NASA/TP-2005-212505, 2005. 
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CARES is a post- 
processor to FEA 
Operates at the 
macro scale of the 
material 



Predicted 
component failure 
probability vs: load 
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Approach for Life Prediction & Component Design of Composites 
> Combine CARES, MAC & FEA codes FEAMAC/CARES 


Move CARES from the macroscopic scale^^ 
of the structure to the microscale of the y 
individual RUC material constituents 

• FEAMAC/CARES Capability: 

• Individual constituent and component level 
probability of failure tracked (for failure initiation) 

• Individual & concurrent failure modes 

• Laminate level analysis capability 

• Progressive damage capability/simulation 

> Subcells killed at random failure thresholds 

Debonding/crack path physics at 
constituent level not explicitly included 
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Fiber Interface Matrix 
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Progressive Damage Criterion 

Calculate failure probability, P f , for each material constjtuent of the RUC 

associated with an element integration point 

CARES calculated RUC 

Pf (CARES) RUC 


Random number generated 


([Random} 


of RUC 


Encourages more rapid 
damage propagation than 
failing individual subcells 
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Fail all material 

constituent subcells 
> Kill elastic modulus 
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Random Element Failure vs: Neighbor Influenced Failure 

(Cellular Automaton Enhancement) 

Encourage more abrupt failure and “ crack-like ” damage growth patterns 

Failed element 

Adjacent element 

Adjacent 2 elements 
with highest 

P f (CARES) has 
P f (Random) adjusted 


A cellular automaton is a 
collection of "colored" cells on 
a grid that evolves through 
discrete time steps according to 
a set of rules based on the 
states of neighboring cells 


Rule: When failure of an 
element is encountered, the 
random failure threshold of the 
neighboring elements are 
adjusted to that of the failed 
element. L oad state determines 
whjch elements have highest 
probability of Mum 


Random element failure 
FF 


Loading 

Direction 



f-fei 


Initial damage 
is diffuse and 
resists 
propagation 


Example: 0° Ply 
uniaxial ramp load 
25x25 FEA mesh 


Cellular automaton 
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Adjusted element P f (Random) 
more likely to be lower than 

original P f(Random) and fail 
sooner as load increases - 
enhancing damage propagation 
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0° single ply tensile specimen (Load parallel to fiber axis) 

Shown are two trial executions; one using the automaton adjusted element 

feature and one with the feature inactive 
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0° single ply tensile specimen 

Progression of damage in FE model of a unidirectional ply under longitudinal loading 


Not Adjusted 



Automaton Adjusted 

m 



I I Matrix failure | Adjacent to failed matrix 
Fiber failure | Adjacent to failed fiber 
No failure 
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(a) and (d ) ; early matrix damage 

(b) and (e ) ; progression to substantial 

matrix damage 

(c) and (f ) ; final composite failure 

(fiber failure) 
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90° single ply tensile specimen (Load transverse to fiber axis) 

Shown are two trial executions; one using the automaton adjusted element 



6rc 


Strain 


Gienn Research center at Lewis Field 39 th International Conference & Expo on Advanced Ceramics and Composites 



90° single ply tensile specimen 

Progression of damage in FE model of a unidirectional ply under transverse loading 
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(a) and (c ) ; early matrix damage 

(b) and (d ) ; final composite failure 

(matrix failure) 
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Example: SiC/RBSN Laminated Composite in On-Axis & Off-Axis Loading 

• Tested by Bhatt & Phillips (1990) 

> displays key mechanisms/features for model material 

• SCS-6 fiber/Reaction Bonded Silicon Nitride matrix composite 
examined in detail by NASA > several papers published 

• Laminated CMCs of interest to industry and less complex than 
woven composites 

> failure modes are not conflicted with complex fiber architecture 

• [0] & [0/90] laminates display nonlinearity due to matrix failure, 
followed by fiber failure. 

• Remaining ply orientations display sudden brittle failure. 


Bhatt, R.T., and Phillips, R.E.: “Laminate Behavior for SiC Fiber-Reiinforced Reaction-Bonded 
Silicon Nitride Matrix Composites.” J. of Comp. Tech. & Res. V. 12, No. 1, Spring 1990, pp. 13-23. 
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Stress (Mpa) 


SiC/RBSN 


700 r— 


Experimental Results 
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SiC/RBSN Example Procedure & Setup 

S 4 Sheiidements Stochastic strength analysis: 

Fixed-displacement ramp load (from individual trials / simulations / realizations) 


1) Cool down 
from stress-free 
temperature of 
550° to room 
temperature 23° 

2) apply fixed- 
displacement 
ramp load 

i 



Loading 

Direction 



❖ 


Use CARES Unit Sphere failure criterion 

assume Isotropic material constituent strength 
- for simplicity and initial testing 


*'?'»*./ Residual stresses I 
l in constituents / 


MAC/GMC RUC 


Weibull parameters correlated 
to experimental results for 0° 
tensile specimen 

Interface strength made large: 

- Encourage matrix to fail 

% before interface 

(10x20 mesh) 

Interfacial failure modes and sliding resistance not considered 




| Fiber 
| Matrix 
Interface 







Constituent properties of SiC/RBSN 

with anisotropic thermal expansion coefficients 


Constituent 

Modulus, GPa 

Poisson ratio 

Longitudinal 
coefficient of 
thermal 

expansion, cll 
(m/m. °C) 

Transverse 
coefficient of 
thermal 
expansion, aj 
(m/m/°C) 

Fiber 

390 

0.17 

4.1xl0" 6 

1.84x1 O’ 6 

Matrix 

110 

0.22 

2.2xl0" 6 

2.2x1 O' 6 

Interface 

1.8 

0.22 

2.0x1 O' 6 

2.0x1 O' 6 


Assumed Weibull Parameters: 

ct oV = 2875 Mpa • m 3/20 
g oV = 150 Mpa • m 3/5 
a oV = 80 Mpa • m 3/5 
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Stress (Pa) 


Residual matrix stresses after cool-down from temperature 

■ Effect of anisotropic fiber-thermal-expansion-coefficient, a f on RUC 
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“ The fiber is orthotropic, with different coefficients of thermal 
expansion along and perpendicular to the fiber axis.”; 

Saigal, A, Kupperman, D. S., Singh, J. P., Singh, D., and J. Richardson (1993): 
“Thermal Residual Strains and Stresses in Silicon Carbide-Fiber-Reinforced Silicon 
Nitride Composites”. Composited Engineering, Vol. 3, No. 11, pp. 1075-1086. 
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0° Sinaia Plv FEAMAC/CARES Calibrated 
y y to 550 0 C stress-free 

temperature & 24% fiber 
volume fraction , V f 


Approx, exp. 
Curve (30% V f ) 


■ Effect of matrix fragments 

On fiber response not modeled 

alibration (24% V f ) | 


ctual stress-strain curve from: 

Chulya, A., Gyekenyesi, J. P., and Bhatt, R. (1991); 
“Mechanical Behavior of Fiber Reinforced 
SiC/RBSN Ceramic Matrix Composites: Theory 
and Experiment. NASATM 103688; 


O CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 

•Exp. 24% fiber vol. 
•Exp. 30% fiber vol. 
50% Fail Prob 24% 

□ 50% Fail. Prob. 30% 

• ••• FEAMAC/CARES 24% 
FEAMAC/CARES 30% 


0 0.001 0.002 0.003 0.004 

Strain 

(GRC 
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0° Single Ply [ 0 ] 

Non-linear (graceful) failure behavior ^ 

jL i 

Calibrated to 
experimental 
data for 24% 
fiber volume 
fraction 


FEAMAC/CARES analysis was for 
a single ply to speed computation 



Exp. 24% fiber vol. 
Exp. 30% fiber vol. 
50% Fail Prob 24% 
50% Fail. Prob. 30% 
FEAMAC/CARES 24% 
FEAMAC/CARES 30% 


Predicted response for 
30% fiber volume fraction 
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FEAMAC/CARES: 
Mean = 238.2 MPa 
Std. Dev. = 10.28 MPa 

Experimental Data: 

Matrix Frac. Str.: 227 +41 Mpa 


• Strain 
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Brewer, D., “ Computation of Variability in the Average Thermal and Mechanical Properties of a Melt-Infiltrated 
smperature Ceramic Matrix Composites 5, M. Singh, R.J. Kearns, E. Lara-Curzio, R. Naslain, Eds, 2004, pp. 279-284 


Damage progression of 0° tensile specimen - two trials (undeformed plot) 
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Prediction for Ten Trials for 90° Fiber Orientation 

Assuming matrix and interface are isotropic strength materials 

•■•Experiment 

□ CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 



1 -of-1 0 
Individual 
simulation 
trials 


Experimental Data: 
Brittle behavior Frac. Str.: 27 +3 MPa 


Note: very few specimens were tested which means the range of 
uncertainty (the confidence bounds ) for m is large/ 


0.0002 


0.0004 0.0006 

Strain 


0.0008 


0.001 


Gienn Research center at Lewis Field 39 th International Conference & Expo on Advanced Ceramics and Composites 



90° Tensile specimen at final failure for 10 trials - Undeformed plots 
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Prediction for Ten Trials for 10° Fiber Orientation 


Assuming matrix and interface are isotropic strength materials 



1 -of-1 0 
Individual 
simulation 
trials 


^^Experi merit 

O CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 


FEAMAC/CARES 
analysis was for a 
single ply to speed 

computation 

[10] 8 


Brittle behavior 

Experimental Data: 
Frac. Str.: 162 MPa 


Note: very few specimens were tested which means the 
range of uncertainty (the confidence bounds ) for m is large/ 
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10° off-axis tensile specimen; 10 trials at final (matrix) failure; deformed plots 

■ Edges are allowed to freely deform (warp) on cool-down 

■ After cool-down; bottom edge fixed in loading direction when displacement load applied 

■ After cool-down; single node along top edge (middle) fixed in direction perpendicular to displacement direct. 
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Pred ction for Ten Trials for 45° Fiber Orientation 

Assuming matrix and interface are isotropic strength materials 

•Experiment 



FEAMAC/CARES 
analysis was for a 
single ply to speed 
computation 



□ CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 


Experimental Data: 

Frac. Str.: 43 MPa 

Note: very few specimens were tested which means the 
range of uncertainty (the confidence bounds ) for m is large/ 
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Prediction for Ten Trials for [+45 2 /-45 2 ] 8 Fiber Orientation 


Assuming matrix and interface are isotropic strength materials 



Matrix cracks 
approx, 
normal to 
loading 
direction 


FEAMAC/CARES 


analysis was for four 
plys (+45/-45M5/+45) to 
speed computation 


[+45 2 /-45p] 


2J8 


jf Neither graceful or brittle behavior 
Experimental Data: 

Matrix Frac. Str.: 75 ±10 MPa; m « 8.8 
Ultimate Str: 88 ±16 MPa; m« 6.3 


•Experiment 

l~l CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 


Note: very few specimens were tested which means the 
range of uncertainty (the confidence bounds ) for m is large/ 
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For [+45 2 1-45 ] Fiber Orientation; 10 trials at final (matrix) failure; 

deformed plots 


FEAMAC/CARES analysis was for four plys (+45/-45M5/+45) to speed computation 
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Prediction for Ten Trials for [0 2 /90 2 ] s Fiber Orientation 

Assuming matrix and interface are isotropic strength materials 
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FEAMAC/CARES 
analysis was for four 
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computation 
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►Experiment 

□ CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 
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Non-linear (graceful) failure behavior 


Experimental Data: 

Matrix Frac. Str.: 127 ±26 MPa 
Ultimate Str: 294 ±87 MPa 

Note: very few specimens were tested which means the 
range of uncertainty (the confidence bounds ) for m is large/ 
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Prediction for Ten Trials for [0 2 /90 2 ] s Fiber Orientation 
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Intersection of 
simulation trial 
with strain offset 
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Experimental Data: 
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0° Double-Notched Tensile Specimen 

Failure mode showed axial splitting of matrix 
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FEAMAC/CARES analysis was for 
a single ply to speed computation 








0° Double-Notched 
vs: Central-Hole 
Tensile Specimen 
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Time-dependent Failure Example: Static Loading 


(Matrix Damage Accumulation From Slow Crack Growth) 
Strain response for applied static tensile load over time 
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stress applied to a 
0° SiC/RBSN ply 

10 time increments 
per time magnitude 

Slow Crack Growth 
Power Law: 

d a _ ait n 

Weibull Parameters 
m = 7 (Weibull slope) 

g 0 = 106 Mpa • mm 3/7 

Fatigue Parameters 
N = 20 (fatigue slope) 
B = 1 .0E9 MPa 2 * sec 



Conclusions 


• Progressive damage simulation of composite structures 
incorporating probabilistic material strength models is possible with 
the FEAMAC/CARES code 

• The Unit Sphere multiaxial model was used predict the strength 
response of a SiC-RBSN composite for various fiber orientations 
under uniaxial tension 

• Reasonable correlation to matrix cracking strength experimental 
data was achieved assuming the matrix was an isotropic material 
with m » 5, and assuming residual stresses from thermal processing 
were present 

• Brittle behavior vs: non-brittle failure (graceful failure) demonstrated 

• Localized damage modes at stress concentration features shown 
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Extra Material 
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Abstract: 

Reported here is a coupling of two NASA developed codes: 
CARES (Ceramics Analysis and Reliability Evaluation of 
Structures) with the MAC/GMC (Micromechanics Analysis 
Code/ Generalized Method of Cells) composite material 
analysis code. The resulting code is called 
FEAMAC/CARES and is constructed as an Abaqus finite 
element analysis UMAT (user defined material). Here we 
describe the FEAMAC/CARES code and an example 
problem (taken from the open literature) of a laminated CMC 
in off-axis loading is shown. FEAMAC/CARES performs 
stochastic-strength-based damage simulation response of a 
CMC under multiaxial loading using elastic stiffness 
reduction of the failed elements. 
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0° single ply tensile specimen (Load parallel to fiber axis) 
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Trade off: 


Mesh sensitivity vs: 
localization of damage 


O 50% Fail. Prob. 

• • • • 10x20 mesh; 320 steps 

• • 10x20 mesh; 100 steps 
•3x6 mesh; 320 steps 
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Time-Dependent Life Prediction Theory - 
Slow Crack Growth and Cyclic Fatigue Crack Growth Laws 


POWGf Law I - Slow Crack Growth (SCG) 



Combined Power Law & Walker Law: SCG and Cyclic Fatigue 
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Time-Dependent Life Prediction Theory - 

Slow Crack Growth and Cyclic Fatigue Crack Growth Laws 

with discrete time steps 

Modeling individual time steps in the life prediction 
methodology enables simulating transient events such as 
turbine start-up/shut-down or atmospheric re-entry. A 
computationally efficient methodology has been 
developed that can extrapolate the reliability calculation 
for an arbitrary number of Z cycles - where each cycle is 
described by k number of time steps. This conceivably 
allows the coupling of other effects such as stiffiness 
degradation and oxidation effects on the individual time 
steps and this can be accounted for interactively within 
the transient finite element and micromechanics analysis. 

6rc 
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Transient Life Prediction Theory - 

Power Law SCG 



Individual time step: Each time step can have different loading, 
Weibull, and fatigue parameters. Compatibility of failure 
probability is maintained between the individual time steps 


SiC/RBSN Notional Example for SCG 

0° Degree tensile specimen under a static load over time 

❖ Use same 10x20 mesh, RUC, and material properties as previous 
SiC/RBSN off-axis loading example 


Weibull and Slow Crack Growth (SCG) Parameters 


Constituent 

Weibull 
modulus, mv 

Weibull scale 
parameter, <j 0 v, 
MPa • mm 3/mF 

Fatigue 
exponent, Nr 
(Equation (11)) 

Fatigue constant, 
IV, MPa 2 • sec 
(Equation (17)) 

Fiber 

20.0 

2875.0 

100.0 

l.OxlO 10 

Matrix 

7.0 

106.0 

20.0 

l.OxlO 9 

Interface 

7.0 

60.0 

100.0 

l.OxlO 10 


Unit Sphere Multiaxial (Batdorf) Model: 

Puts linear elastic fracture mechanics into Weibull weakest-link theory 


> Incremental failure probability is 
the product of two probabilities: 

A Pf =PyP2 

Pi = Probability of the existence of a crack 
having a critical strength between a c and a c + 

A a c in the incremental volume 

P2 = Probability a crack having a critical strength 

of g c will be oriented in a direction such that it will 
fail under the applied multiaxial stress state 


P2 involves Integration of an equivalent stress a e , 
where a e >a c , over the surface of a unit radius sphere 
(all possible flaw orientations) divided by the total surface 
area of the unit radius sphere 



a e is a function of an assumed crack shape and 
multiaxial fracture criterion 

Mixed-Mode Fracture Criteria: 


> Component failure probability: 





• Normal stress (shear-insensitive cracks) 

• Maximum tensile stress 

• Total coplanar strain energy release rate 

• Noncoplanar (Shetty) 

Flaw Shapes: 
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Griffith crack 
Penny-shaped crack 
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CARES Unit Sphere Multiaxial model 

has crack geometry & mixed-mode fracture criterion 

> Two models for transverse isotropy 


(1) Flaw / Fracture-Plane Orientation Anisotropy 


“Polar Cap” (L) 


Polar cap PDF - 
Probability | 
Density 
Distribution f 
Function J 
with shape 
Exponent, (f) L 


Probability 
of occurence 
(frequency) 



“Equatorial 
Belt” (T) 


All possible 
vectors 
normalto / 
the flaw / 
planes 



r- Equatorial 
j bolt PDF 
with shape 
Exponent, <f) T 


Isotropic: 

Random flaw 
orientation 



Direction normal 
to the crack plane 

dA = sin a da dp 


Unit Sphere 


(2) Strength Orientation Anisotropy 
(7, c or K/c varies with orientation 



Nemeth, N.N. (2014): Unit-sphere multiaxial stochastic- 
strength model applied to a composite material. J. 
Comp. Mat., Vol. 48(27) Nov. 2014, pp. 3395-3424 
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Anisotropic Unit Sphere model defined in a 
material coordinate system reference frame 


All possible vectors normal 
to the flaw plane of the fiber 
or tow or unidirectional 
lamina 




Width of belt defines 
extent of misalignment 


x 



All possible flaw plane / 
fracture plane 
orientations of the fiber 
or tow 



►a. 


Similar to Puck’s composite 
failure criterion except in a 
probabilistic framework 

When represented on a unit radius sphere of 
all possible flaw plane orientations the flaw 
plane normals define an equatorial belt 
distribution 
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Multiaxial Performance l biaxial response predicted from a 
MAC/GMC RUC for combined longitudinal (L) and transverse (T) 
loading on a unidirectional PMC vs: FEA. 

50% probability of failure envelope . 

GMC RUC Used 









Unit Sphere parameters 
adjusted so GMC results 
matched FEA results for 
uniaxial tension and 
compression. 

Intermediate points are 
predictions 
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Differences in RUC stress fields 
for a transverse strain: 

GMC RUC 


1 
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HF ■ GMC RUC 
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45° off-axis tensile specimen; 10 trials at final failure deformed plots 

■ Edges are allowed to freely deform (warp) on cool-down 

■ After cool-down; bottom edge fixed in loading direction when displacement load applied 

■ After cool-down; single node along top edge (middle) fixed in direction perpendicular to displacement direct. 
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Prediction for Ten Trials for [0 2 /90 2 ] s Fiber Orientation 

Closer view 


CL 130 



Experiment 

□ CARES calculated 
50% matrix failure 
probability prior to 
any damage initiation 
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0.00075 0.00085 0.00095 0.00105 0.00115 

Strain 
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For [0 2 /90 ] s fiber orientation; four trials with deformed plots 


Progression from matrix failure to final fiber failure 


FEAMAC/CARES analysis was for four 
plys (0/90/90/0) to speed computation 
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Path Forward 

• Continue demonstrate/benchmark capability on CMCs (using available literature) 

> For uniaxial & multiaxial failure response 

(orientation, lamination, stress concentration, flexural) 

- Fast-fracture 

- Time & cycle dependent 

- more detailed micromechanical models of failure modes 

• Develop / incorporate enviromechanical degradation models 

• Investigate applicability to predict EBC damage progression 

• Develop / incorporate anisotropic elastic modulus degradation based on 
CARES critical fracture angle probability density distribution 

• Improve software efficiency (memory, speed, multiprocessing) 

• Demonstrate this capability on component/structure 
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Approach For Life Prediction & Component Design Of Composites 

> Combine CARES, MAC & FEA codes where 

Micromechanics provides the link between structures & materials 


CARES: monlithic ceramics 


• Probabilistic strength 

• Mechanistic-based multiaxial failure model 

• Efficient life prediction algorithm 

• Isotropic and transverse isotropy 

MAC/GMC: composites analysis 



Micromechanics 
Accurate RUC stress fields 
Flexibility in RUC designs 
Progressive damage capability 
Computationally efficient 
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Move CARES from the macroscopic scale of 
the structure to the microscale of the 
individual material constituents & RUC 
with the FEA-MAC micromechanics code 




FEAMAC/CARES 


(Abaqus UMAT) 


Structural-Scale FEA 

Element/Integration 

Point 


Micromechanics Analysis 

(MAC/GMC) 


(CARES) 

Subroutine 
Reliability analysis at 
the RUC level 


4 a- — -.a A 

Fiber Interface Matrix 
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